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The  Low  Background  Infrared  Calibra- 
tion Facility  (LBIR)  at  the  National  In- 
stitute of  Standards  and  Technology 
has  been  in  operation  for  calibration 
measurements  of  the  radiant  power 
emitted  from  infrared  radiation  (IR) 
sources,  such  as  cryogenic  blackbodies, 
for  more  than  2 years.  The  IR  sources 
are  sent  to  NIST  by  customers  from  in- 
dustry, government,  and  university  labo- 
ratories. An  absolute  cryogenic 
radiometer  is  used  as  the  standard  de- 
tector to  measure  the  total  radiant 
power  at  its  aperture.  The  low  back- 
ground is  provided  by  a closed  cycle 
helium  refrigeration  system  that  main- 
tains the  inner  parts  of  the  calibration 


chamber  at  20  K.  The  radiance  temper- 
ature of  the  blackbody  is  deduced  from 
the  measured  power  and  compared 
with  the  blackbody  temperature  sensor 
data.  The  calibration  procedures  and 
data  analysis  are  illustrated  using  the 
measurements  of  a typical  blackbody. 
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1.  Introduction 

The  national  need  to  establish  a primary  stan- 
dard reference  for  calibrating  infrared  sensors  led 
to  the  establishment,  in  1989,  of  the  Low  Back- 
ground Infrared  Calibration  Facility  (LBIR)  at  the 
National  Institute  of  Standards  and  Technology. 
This  facility  has  been  serving  a variety  of  users  from 
industry  and  government  laboratories  by  providing 
calibrations  of  cryogenic  blackbody  sources.  These 
sources  are  used  by  the  customer  for  calibrations  of 
their  infrared  sensor  systems.  An  absolute  cryo- 
genic radiometer  (ACR)  of  the  electrical  substitu- 
tion type  has  been  developed  as  the  standard 
reference  detector  for  the  LBIR  calibrations.  The 
ACR  measures  the  total  radiant  power  incident  on 
its  precision  aperture.  Its  operating  range  for  radi- 
ant power  measurements  is  20  nW  to  100  (xW  with 


an  expanded  uncertainty  of  ±1%  (95%  confidence 
level).  The  physical  description  of  the  ACR  and  its 
characterization  as  an  absolute  detector  for  radiant 
power  measurements  have  been  described  in  Ref. 
[1].  The  LBIR  chamber  is  shown  in  Fig.  1. 

The  cylindrical  vacuum  tank  housing  the  ACR 
detector  and  the  customer  blackbody  source  is 
made  of  stainless  steel  and  is  152  cm  long  and  60  cm 
in  diameter.  Two  shields  made  of  copper  (shown  in 
the  cutaway  portion  of  Fig.  1)  separated  by  a 2.54 
cm  gap  are  cooled  by  a helium  refrigerator  which 
circulates  cooled  helium  gas  (15  K)  through  copper 
lines  vacuum  brazed  to  these  shields.  The  inner 
shield  operates  at  20  K and  provides  the  low  ther- 
mal radiation  background  inside  the  chamber.  The 
inner  surface  of  the  inner  shield  has  a highly 
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Fig.  1.  LBIR  chamber  with  partial  cutaway  showing  the  major  features  of  the  apparatus. 


absorbent  black  coating  of  IITRI  MH-2200  paint' 
to  reduce  scattered  radiation.  The  ACR  can  be 
placed  in  any  of  the  three  ports  of  the  LBIR  cham- 
ber, which  are  located  at  distances  of  approxi- 
mately 30  cm,  65  cm,  and  100  cm  from  the 
blackbody  source  under  test.  The  placement  of  the 
ACR  can  be  chosen  so  that  the  radiant  power  at  its 
aperture  is  in  the  ACR  measurement  range.  A 
complete  description  of  the  facility  can  be  found  in 
a previous  publication  [3]. 

The  general  requirements  for  blackbody  calibra- 
tions are  described  in  Sec.  2 of  this  article.  The 
procedures  for  handling  the  cryogenics  of  the  cali- 
bration setup  and  the  procedures  adopted  for  mea- 
surements and  analysis  are  described  in  Sec.  3.  A 
recent  blackbody  calibration  illustrates  the  proce- 
dures. The  blackbody  in  question  has  two  sensors 
and  is  referred  to  as  BB  in  this  paper. 

2.  Calibration  Requirements 

Calibration  requirements  fall  into  two 
categories:  1)  needs  established  by  the  customer 
and  2)  needs  of  the  NIST  LBIR  facility.  A typical 


' Certain  commercial  equipment,  instruments,  or  materials  are 
identified  in  this  paper  to  specify  adequately  the  experimental 
procedure.  Such  identification  does  not  imply  recommendation 
or  endorsement  by  the  National  Institute  of  Standards  and 
Technology,  nor  does  it  imply  that  the  materials  or  equipment 
identified  are  necessarily  the  best  available  for  the  purpose. 


customer  requirement  for  the  BB  is  as  follows:  1) 
measure  the  radiant  power  at  the  ACR  aperture 
for  the  blackbody  temperature  range  of  200  K to 
400  K in  increments  of  25  K for  the  source  aper- 
ture diameter  of  650  |xm  and  2)  calibrate  the  tem- 
perature sensors  of  the  source  by  comparing  with 
the  deduced  radiance  temperature.  NIST  issues 
guidelines  to  the  customer  for  preparing  the  black- 
body source  to  conform  to  the  vacuum  and  allow- 
able contaminant  level  specifications  of  the  LBIR 
chamber.  The  NIST  requirement  specifies  that  the 
partial  pressure  of  hydrocarbons  (beyond  approxi- 
mately 45  u)  should  not  exceed  1.33  x 10"*  Pa  for  a gj 
total  base  pressure  of  1.33x10"'’  Pa.  Other  re- 
quirements are  described  in  detail  in  Ref.  [3]. 

NIST  also  supplies  the  customer  with  a mounting  jl 
plate  and  guidelines  for  installing  the  plate  on  the  I 
blackbody  source  [3]. 

3.  Calibration  Procedures 

3.1  Cryogenic  Apparatus  and  Preparation  for 
Measurements 

The  blackbody  is  integrated  into  the  LBIR 
chamber  after  testing  its  vacuum  integrity.  Cryo- 
conditioning  of  the  LBIR  chamber  follows  with  the 
inner  shield  of  the  chamber  being  cooled  to  ap- 
proximately 20  K and  the  outer  shield  being  cooled 
to  approximately  40  K.  The  blackbody  mounting 
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shelf  and  the  main  chamber  have  separate  cold  he- 
lium feeds  with  flow  control  valves  to  allow  diver- 
sion of  additional  coolant  to  the  blackbody,  if 
needed.  The  cryoshields  have  several  silicon  diode 
thermometers  attached  at  various  points.  These 
thermometers  are  monitored  during  calibration  to 
ensure  stability  of  the  IR  environment. 

In  order  to  provide  a steady  background  flux  for 
the  ACR,  a copper  disk  separates  the  blackbody 
chamber  from  the  ACR  chamber.  This  disk,  called 
the  isothermal  plate,  is  moderately  thermally  insu- 
lated from  the  inner  shield.  The  isothermal  plate  is 
fitted  with  an  aperture  that  blocks  the  background 
radiation  from  the  blackbody  chamber  without  lim- 
iting the  field  of  view  of  the  ACR  aperture.  Com- 
plete isolation  of  the  ACR  chamber  is  achieved,  if 
necessary,  by  a shutter  in  front  of  the  isothermal 
plate  that  can  be  operated  by  remote  control.  The 
shutter  is  thermally  anchored  to  the  isothermal 
plate  with  a flexible  copper  braid.  The  temperature 
of  the  isothermal  plate  and  shutter  are  monitored 
by  three  precision  silicon  diode  thermometers.  The 
isothermal  plate  is  controlled  to  0.05  K by  using  an 
integral  heating  element  and  the  temperature  is 
held  at  22  K. 

3.2  Measurements 

Figure  2 shows  the  calibration  setup  inside  the 
LBIR  chamber.  Blackbody  sensor  data  and  radiant 
power  data  are  collected  at  each  blackbody  tem- 
perature and  aperture  setting  specified  by  the  cus- 
tomer. Data  are  collected  and  recorded  at 
approximately  one  second  intervals  for  a period  of 
three  minutes  at  each  setting.  These  data  form  the 
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Fig.  2.  Blackbody  calibration  setup  inside  the  LBIR  chamber. 


basis  for  estimating  the  repeatability  of  the  instru- 
mentation. Cryogenic  preparation  of  the  ACR  and 
the  complete  calibration  sequence  are  repeated  on 
three  different  days  to  test  the  reproducibility  of 
the  entire  measurement  system. 

3.2.1  Blackbody  Temperature  Sensors  The 
blackbody  temperature  is  generally  measured  by 
platinum  resistance  thermometers  (PRTs) 
mounted  on  the  blackbody  core.  Voltages  across 
the  PRTs  measured  by  a multimeter,  and  currents 
supplied  by  a constant  current  source  are  con- 
verted to  resistance  values  using  Ohm’s  law.  The 
resistances  are  converted  to  temperature  values  if 
the  customer  provides  data  on  the  calibration  of 
the  PRT  sensors.  Nominal  settings  on  the  external 
temperature  control  electronics  of  the  BB  and 
measured  temperatures  for  the  BB  sensors  are 
shown  in  columns  1,  2,  and  4 of  Table  1.  Experi- 
mental standard  deviations  from  the  1 s repetitions 
at  each  setting  on  each  day  are  shown  in  columns  3 
and  5 of  the  same  table. 

3.2.2  Radiant  Power  The  radiant  power  at 
the  ACR  aperture  is  measured  as  follows:  The  ra- 
diation from  the  blackbody  is  blocked  from  the 
ACR  receiver  by  using  the  shutoff  position  on  the 
blackbody  aperture  wheel  (if  available).  The  elec- 
trical heater  power  to  the  ACR  receiver  is  set  to  a 
value  higher  (preferably  20%  to  30%)  than  the  ex- 
pected radiant  power  from  the  blackbody  so  that 
the  temperature  controller  servo  (ac  bridge)  can 
maintain  a constant  temperature  at  the  ACR  cone 
by  varying  the  amount  of  electrical  power.  Prelimi- 
nary estimates  of  radiant  power  at  the  ACR  aper- 
ture are  made  using  the  Stefan-Boltzmann  law  and 
the  approximate  distance  between  the  ACR  aper- 
ture and  the  blackbody  aperture.  After  the  ac 
bridge  is  balanced,  the  ACR  is  allowed  to  receive 
the  blackbody  radiation  by  turning  the  aperture 
wheel  to  a chosen  aperture  position.  If  the  cus- 
tomer blackbody  is  equipped  with  a fixed  aperture, 
the  152  mm  (6  in)  shutter  in  front  of  the  isothermal 
plate  is  used  for  blocking  the  blackbody  radiation. 
The  difference  between  the  initial  electrical  power 
setting  with  the  shutter  closed  and  the  final  electri- 
cal power  setting  with  the  shutter  open  is  the  radi- 
ant power  at  the  ACR  aperture.  Data  is  taken 
continuously  at  one  second  intervals  for  at  least  3 
min  with  the  shutter  open,  in  order  to  determine 
the  repeatability  of  the  instrumentation. 

3.23  Geometry  The  measurement  of  the  ra- 
dius of  the  blackbody  aperture,  n (see  Fig.  2),  cor- 
rected for  thermal  contraction  at  20  K,  is  supplied 
by  the  customer.  The  radius  of  the  ACR  aperture, 
ri,  has  been  measured  by  the  Precision  Engineering 
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Table  1.  Experimental  data 


Nominal 

(K) 

Sensor  1 
(K) 

Std.  dev. 

(K) 

Sensor  2 
(K) 

Std.  dev. 

(K) 

Measured 

power(nW) 

Average 
power  (nW) 

Std.  dev.  of 
average  (nW) 

200 

199.897 

0.002 

199.281 

0.002 

71.9 

71.5 

0.5  (0.6%) 

199.874 

0.005 

199.332 

0.005 

72.0 

199.947 

0.002 

199.328 

0.003 

70.6 

225 

224.780 

0,003 

224.266 

0.004 

113.5 

113.3 

0.1  (0.1%) 

224.728 

0.003 

224.276 

0.006 

113.4 

224.823 

0.004 

224.306 

0.004 

113.0 

250 

249.662 

0.004 

249.268 

0.004 

171.3 

172.0 

0.4  (0.2%) 

249.639 

0.007 

249.288 

0.008 

172.6 

249.733 

0.005 

249.322 

0.004 

172.2 

275 

274.685 

0.003 

274.408 

0.002 

252.6 

252.5 

0.2  (0.1%) 

274.644 

0.006 

274.394 

0.007 

252.2 

274.719 

0.006 

274.443 

0.004 

252.8 

300 

299.531 

0.007 

299.329 

0.007 

358.4 

358.0 

0.4  (0.1%) 

299.521 

0.003 

299.346 

0.004 

357.3 

299.588 

0.005 

299.389 

0.004 

358.4 

325 

324.412 

0.005 

324.346 

0.006 

497.6 

494.6 

1.4  (0.3%) 

324.457 

0.005 

324.375 

0.005 

492.8 

324.446 

0.003 

324.379 

0.004 

493.4 

350 

349.306 

0.003 

349.378 

0.003 

665.8 

664.8 

0.8  (0.1%) 

349.337 

0.008 

349.395 

0.011 

663.7 

349.445 

0.006 

349.479 

0.006 

664.8 

375 

374.047 

0.005 

374.284 

0.005 

882.2 

878.3 

2.0  (0.2%) 

374.069 

0.003 

374.292 

0.003 

875.7 

374.171 

0.006 

374.377 

0.007 

877.0 

400 

399.026 

0.005 

399.362 

0.005 

1138.0 

1137.1 

1.3  (0.1%) 

399.056 

0.004 

399.379 

0.004 

1134.6 

399.131 

0.004 

399.444 

0.003 

1138.8 

Division  at  NIST  at  room  temperature  in  three  dif- 
ferent orientations.  The  average  radius  is  corrected 
for  thermal  contraction  as  explained  in  Ref.  [1]. 

The  distance,  S,  between  the  blackbody  aperture 
and  the  ACR  limiting  (precision)  aperture  is  mea- 
sured in  two  parts.  The  user  measures  the  distance 
between  the  blackbody  aperture  and  the  front  sur- 
face of  a reference  tab  on  the  blackbody  mounting 
plate  as  required  by  NIST  [3],  and  NIST  personnel 
measure  the  distance  from  the  front  surface  of  the 
tab  to  the  ACR  aperture.  The  chamber  measure- 
ment is  assisted  by  a Kaman  proximity  sensor  [4] 
located  internally  to  determine  the  final  location  of 
the  reference  tab  to  the  ACR  aperture  before  the 
chamber  is  evacuated.  The  measurements  are 
made  at  ambient  temperature  and  corrected  for 
thermal  contraction  due  to  cooling  of  the  chamber 
and  the  ACR  to  20  K and  2 K,  respectively.  The 
radius  of  the  aperture  in  the  isothermal  plate,  rs. 


shown  in  Fig.  2,  is  chosen  to  limit  the  background 
flux  from  the  blackbody  front  surface.  However,  it 
will  allow  the  full  cone  of  light  from  the  blackbody 
aperture  to  reach  the  ACR  aperture.  The  geomet- 
rical measurements  for  the  BB  are  given  in  Table 
2. 

3.2.4  Radiance  Temperature  The  following 
equation  deduced  from  the  Stefan-Boltzmann  law 
is  used  to  convert  the  radiant  power  to  radiance 
temperature,  T: 


T = \.F\A\  aM. 

where 

Fi  = (1/2)[2-  [2^-4a:^ , 

x=|,y  =^,z  = l + (l+jc')y^ 
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Table  2.  Geometric  data  for  BB  calibration 

Meas. 

Std.  dev. 

value 

mean 

1. 

Radius  {ri)  of  the  ACR  aperture  at  2.2  K 

1.4971  cm 

5 X 10”®  cm 

2. 

Distance  between  the  apertures  (A) 

30.77  cm 

0.042  cm 

3. 

BB  aperture  radius  (/-j): 

Measured  at 

Deduced 

Nominal 

ambient 

20  K 

Std.  dev.  mean 

mm 

mm 

mm 

mm 

0.325 

0.3251 

0.3244 

0.0013 

andTi  = 'nTi^  is  the  area  of  the  blackbody  aperture; 
P is  the  radiant  power  and  ctm  is  the  Stefan-Boltz- 
mann  constant,  and  equal  to  5.67051  x 10“**  W m“^ 
K”"*.  The  expression  for  the  configuration  factor, 
F:,  given  above  is  taken  from  Ref.  [5]  and  to  a first 
order  approximation  reduces  to  The  quan- 

tities ri,  vi,  and  S have  already  been  defined,  and 
their  measured  values  at  room  temperature  are  ad- 
justed to  cryogenic  temperatures  using  the  stan- 
dard reference  data  for  contraction  of  materials 
[6]. 

3.3  Analysis 

3.3.1  Diffraction  Correction  It  has  been  es- 
tablished in  the  literature  [7,  8]  that  the  radiation 
from  the  blackbody  incident  at  the  ACR  aperture 
will  not  be  solely  determined  according  to  geomet- 
rical optics  because  of  diffraction  effects  at  both 
limiting  and  nonlimiting  apertures  in  the  beam 
path.  Therefore,  diffraction  losses  at  each  one  of 
the  apertures  in  the  beam  path  are  estimated  by 
using  the  procedures  published  in  Refs.  [7]  and  [8]. 
The  measured  radiant  power  values  and  the  de- 
duced radiance  temperatures  are  corrected  for  the 
total  diffraction  loss  which  is  a product  of  the  frac- 
tional losses  at  each  aperture  in  the  beam  path. 
The  calculated  correction  factors  (AR/R)  given  as  a 
percentage  of  the  measured  power  values  for  the 
BB  are  listed  in  Table  3 for  various  temperatures. 
The  corresponding  correction  factor  (AT/T)  for 
the  radiance  temperature  is  approximately  1/4  of 
the  correction  (AP/F)  for  the  radiant  power. 


Table  3.  Calculated  diffraction  correction  (AP)  given  as  a per- 
centage of  the  measured  radiant  power  (P).  Systematic  uncer- 
tainty (Type  B)  in  the  correction  due  to  approximations  in  the 
calculations  is  ±10%  of  AP/P. 


Nominal 
blackbody 
temp.  ( K ) 

Diff.  corr. 

A P/P 

Uncer.  (1  it) 

200 

1.3% 

±0.13% 

225 

1.2% 

±0.12% 

250 

1.1% 

±0.11% 

275 

1.0% 

±0.1% 

300 

0.9% 

± 0.09% 

325 

0.8% 

±0.08% 

350 

0.8% 

±0.08% 

375 

0.7% 

±0.07% 

400 

0.7% 

±0.07% 

3.3.2  Measurement  Uncertainties  The  uncer- 
tainties in  the  measurements  are  analyzed  accord- 
ing to  guidelines  under  development  by  the 
International  Organization  for  Standardization 
(ISO)  [9].  Uncertainties  are  caused  by  either 
systematic  or  random  effects.  Uncertainties  from 
systematic  effects  (here  after  called  systematic  un- 
certainties for  simplicity)  are  of  the  following  types 
for  this  experiment:  1)  uncertainties  in  theoretical 
corrections  applied  to  measured  data  such  as  the 
uncertainties  associated  with  the  diffraction  correc- 
tions; 2)  uncertainties  in  the  measured  values  of 
parameters  that  remain  constant  during  the 
calibration.  Examples  of  such  uncertainties  are  the 
parameters  that  characterize  the  ACR  response  [1] 
which  do  not  change  during  the  ACR  power 
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measurements.  These  uncertainties  are  Type  B un- 
certainties according  to  Ref.  [9]  Subclause  3.3.3. 
Uncertainties  from  random  effects  (hereafter 
called  random  uncertainties  for  simplicity)  arise 
from  the  statistical  variation  in  the  measurements 
during  the  calibration  experiment.  This  type  of 
variation  is  inherent  in  the  blackbody  PRT  sensor 
measurements  of  temperature  and  the  ACR  mea- 
surements of  radiant  power  during  the  calibration 
experiment.  These  are  Type  A uncertainties  [9]. 

Uncertainty  in  Radiant  Power 

The  uncertainty  in  the  ACR  measurement  of  the 
radiant  power  is  analyzed  as  follows:  For  each 
nominal  temperature,  radiant  power  data  (F,)  are 
collected  every  1-1/2  s for  approximately  3 min, 
leadiiig  to  approximately  120  data  points.  An  aver- 
age, P,  and  standard  deviation,  5rpt,  are  calculated 
from  this  data.  The  standard  deviation,  5rpt,  is  a 
measure  of  the  repeatability  of  the  measurements 
in  each  set,  and  in  general,  the  ACR  data  show  very 
good  repeatability.  The  measurement  process  is  re- 
peated three  different  times  to  estimate  a standard 
deviation,  5rpr,  for  reproducibility.  The  average  over 
the  three  runs  is  calculated  as: 

P = \ IPj 

a standard  deviation  (5t)  is  calculated  from 

s;='Z(h-pyn.  (4) 

7=1 

The  standard  deviation  (^t)  is  a measure  of  both  re- 
peatability and  reproducibility  of  measurements 
[10]  according  to 

s\~s^r+s]ftl\2Q . (5) 

Although  more  repetitions  of  the  entire  measure- 
ment process  would  be  desirable  for  better  statis- 
tics, time  constraints  limit  _the  number  of 
repetitions.  The  radiant  powers,  Pj,  and  the  average 
powers,  P,  for  the  BB  are  given  in  columns  6 and  7, 
respectively,  of  Table  1.  The  standard  deviations  of 
the  means,  Jt/\/3,  are  given  in  column  8. 

Two  sources  of  systematic  uncertainty  are  evalu- 
ated as  Type  B following  Ref.  [9].  The  standard  un- 
certainty (i.e.,  1 O' level)  in  the  diffraction  correction 
to  the  measured  radiant  power  is  based  on  diffrac- 
tion theory  used  to  calculate  the  correction.  This 
standard  uncertainty  is  estimated  to  be  ± 10%  [7,  8] 
of  the  calculated  correction  based  on  scientific 


judgment  (Sec.  4.6  Ref.  [11]),  as  shown  in  column  3 
of  Table  3 for  the  BB.  The  standard  uncertainty  in 
the  characterization  of  the  radiometer  as  given  in 
Ref.  [1]  is  0.12%  of  the  measured  radiant  power. 
The  standard  uncertainty  due  to  systematic  effects 
is  then  the  square  root  of  the  sum  of  the  squares  of 
the  two  components. 

The  random  and  systematic  uncertainties  are 
listed  separately  in  the  calibration  report.  The  com- 
bined standard  uncertainty  (uc)  is  reported  as  the 
square  root  of  the  sum  of  the  squares  of  all  of  the 
components.  However,  it  should  be  noted  that 
other  methods  of  combining  components  of  uncer- 
tainty from  systematic  and  random  effects  are  in 
practice  for  calculating  the  total  uncertainty  [12]. 
Therefore,  all  components  of  uncertainty  are  listed 
to  give  the  customer  the  choice  of  using  an  alterna- 
tive method. 

The  radiant  power  data  from  Table  1 after  cor- 
rection for  diffraction  loss  are  shown  in  column  6 of 
Table  4.  The  uncertainty  components  associated 
with  the  average  power  are  shown  in  columns  7 and 
8,  and  the  combined  standard  uncertainty,  Uc,  is 
given  in  column  9. 

Uncertainty  in  the  Blackbody  Sensor  Data 

The  Type  A (random)  uncertainty  in  the  black- 
body PRT  sensor  data  is  analyzed  as  follows:  The 
PRT  sensor  data  are  collected  in  conjunction  with 
the  ACR  power  measurements.  In  general,  the 
standard  deviations  for  the  PRT  sensor  data  are  ho- 
mogeneous and  negligibly  small.  Repeated  sets  of 
data  are  averaged  to  obtain  the  mean  and  the  corre- 
sponding standard  deviation  analogously  to  the 
method  discussed  earlier  for  the  ACR  power  mea- 
surements. The  mean  temperature  values  of  the  BB 
with  associated  standard  deviations  are  shown  in 
columns  2 through  5 of  Table  4. 

Radiance  Temperature  and  Calibration  Uncertainty 

The  final  radiance  temperatures  for  the  calibra- 
tion report  are  obtained  by  the  least-squares  analy- 
sis of  temperatures,  T/y  (/  = 1,  ; =1,  . . .,3), 

deduced  from  Eq.  (1),  as  a function  of  blackbody 
sensor  readings,  X,j  (/  = 1, . . . , n;  y = 1, . . . , 3),  for 
each  aperture  setting. 

In  order  to  evaluate  a confidence  band  for  the 
variability  of  the  calibration  curve,  the  following 
statistical  procedure  is  adopted.  The  calibration 
equation  is  assumed  to  follow  the  model 

T,j  =flo  + fli  Xij  +fl2  X}j  +a3  Xfj  + ....-F 

OkX^+Cij,  (6) 
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Table  4.  Average  blackbody  sensor  data  and  radiant  jxjwer  corrected  for  diffraction  effects  at  the  ACR  aperture  for  the 
BB 


Nominal 

(K) 

Sensor  1 
data 
(K) 

Std.  dev. 

(K) 

Sensor  2 
data 

(K) 

Std.  dev. 

(K) 

Average 

power 

nW 

Syst. 

uncer. 

nW 

Random 

uncer. 

nW 

Comb. 

uncer. 

Uc  (1  <r) 
nW 

200 

199.906 

0.003 

199.313 

0.003 

72.4 

0.1 

0.5 

0.5  (0.6%) 

225 

224.777 

0.003 

224.282 

0.004 

114.6 

0.2 

0.1 

0.2  (0.2%) 

250 

249.678 

0.005 

249.293 

0.005 

173.9 

0.3 

0.4 

0.5  (0.3%) 

275 

274.683 

0.005 

274.415 

0.005 

255.0 

0.4 

0.2 

0.5  (0.2%) 

300 

299.547 

0.005 

299.355 

0.005 

361.2 

0.5 

0.4 

0.6  (0.2%) 

325 

324.429 

0.005 

324.367 

0.005 

498.6 

0.7 

1.4 

1.6  (0.3%) 

350 

349.362 

0.006 

349.429 

0.007 

669.8 

1.0 

0.8 

1.3  (0.2%) 

375 

374.120 

0.005 

374.331 

0.005 

884.5 

1.2 

2.0 

2.3  (0.3%) 

400 

399.071 

0.005 

399.395 

0.004 

1144.6 

1.6 

1.3 

2.1  (0.2%) 

where  Oo,. . . , a*  are  to  be  estimated,  and  €,j,  the 
random  errors  associated  with  the  measurements 
Tu.  are  assumed  to  be  independent  with  heteroge- 
neous variances  aij.  The  aij  are  estimated  by  s,y 
from  approximately  120  data  points  for  each  run. 
Weighted  least  squares  [13]  accounts  for  the  het- 
erogeneity of  variances  where  the  weights,  w,  , are 
calculated  from  the  empirical  variances  by 

Wij  = l/s,j  . (7) 

The  standard  deviations  associated  with  the  Tij 
are  an  order  of  magnitude  larger  than  the  standard 
deviations  associated  with  the  Xjj  data,  and,  there- 
fore, random  uncertainties  associated  with  the 
measurements  Xij  are  assumed  to  be  negligible  for 
the  purpose  of  the  least  squares  analysis.  The  de- 
gree, k,  for  the  polynomial  in  Eq.  (6)  is  determined 
by  a goodness-of-fit  test  [14]  which  compares  the 
agreement  among  groups  of  three  runs  with  the 
overall  fit  to  the  data.  The  lowest  degree  polyno- 
mial which  satisfies  the  goodness-of-fit  criterion  is 
taken  as  the  calibration  curve. 

Given  a future  blackbody  sensor  setting,  Xh,  its 
calibrated  radiance  temperature  value  is  given  by 

Th  = flo  + Xh  +fl2  Xh  +fl3  Xh  + . . . . Xh  , (8) 

where  ah, , a*  are  least-squares  estimates  from 
calibration  data.  The  random  component  of  uncer- 


tainty associated  with  the  predicted  value  is  com- 
puted as 

f/random  = [(A:  + l)E(95;A:-Fl, 

n-k-\)]'^s{Th)  . (9) 

The  constant  F(95;  ^-l-l,Af-A:-l)is  the  upper  95 
percentile  of  Snedecor’s  F-distribution  with  k + \ 
degrees  of  freedom  in  the  numerator  and  n -k-\ 
degrees  of  freedom  in  the  denominator,  and  s (Th) 
is  the  standard  deviation  of  the  predicted  value,  Th. 
This  uncertainty,  based  on  the  Working-Hoetelling 
confidence  bands  for  the  calibration  curve,  is  valid 
for  all  future  applications  of  the  calibration  curve 
as  long  as  the  model  holds  [15].  The  details  of  the 
statistical  analysis  can  be  found  in  Ref.  [15,  16].  A 
software  package,  OMNITAB  [17],  is  used  for  the 
statistical  computations. 

Data  for  the  least-squares  analysis  of  the  BB  are 
shown  in  Table  5.  Temperature  as  measured  by 
PRT  sensor  1 is  used  as  the  independent  variable, 
X.  The  analysis  confirms  that  a linear  function 
{k  = 1)  is  sufficient  for  describing  the  data,  and  the 
following  equation 

Th=ao+a[(Xh)  (10) 

gives  the  predicted  radiance  temperature  (Th)  for  a 
blackbody  temperature  (Xh)  measured  by  PRT 
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Table  5.  Data  for  the  least-squares  analysis  of  the  BB 


Nominal 

(K) 

Sensor  1 
(K) 

X 

Radiance 

temperature 

(K) 

T 

Standard 

deviation 

(K) 

s 

Weight 

W 

200 

199.897 

201.427 

3.500 

0.0798 

199.874“ 

201.539“ 

1.420 

0.4957 

199.947 

200.546 

2.130 

0.2203 

225 

224.780 

225.736 

2.990 

0.1116 

224.728 

225.666 

1.384 

0.5224 

224.823 

225.506 

1.975 

0.2564 

250 

249.662 

250.135 

2.100 

0.2278 

249.639 

250.593 

1.166 

0.736 

249.733 

250.463 

1.865 

0.2875 

275 

274.6849 

275.575 

1.102 

0.8239 

274.644 

275.466 

1.319 

0.5751 

274.719 

275.605 

1.101 

0.8244 

300 

299.5309 

300.683 

0.715 

1.9545 

299.5211 

300.462 

1.135 

0.7760 

299.588 

300.683 

0.887 

1.2702 

325 

324.412“ 

326.342“ 

0.830 

1.4529 

324.457 

325.600 

0.603 

2.7502 

324.446 

325.658 

1.076 

0.8641 

350 

349.306 

350.931 

0.644 

2.4084 

349.337 

350.657 

0.997 

1.0056 

349.445 

350.807 

0.869 

1.3233 

375 

374.047“ 

376.500“ 

0.660 

2.2930 

374.069 

375.766 

0.816 

1.5003 

374.171 

375.915 

1.467 

0.4647 

400 

399.026 

401.171 

0.872 

1.3164 

399.056 

400.866 

1.266 

0.6242 

399.131 

401.240 

1.217 

0.6748 

“ Data  with  standardized  residuals  from  the  linear  fit  exceeding  2.5.  These  outliers  are  not  included  in 
the  final  analysis. 


sensor  1.  The  estimated  coefficients  and  associated 
standard  deviations  are: 

flo=  -1.096;  5fflo)  = 0.244 
fl'f=  1.007;  s(fl;)  = 0.001  . 

The  percentage  uncertainty  in  the  deduced  radi- 
ance temperature,  8TIT,  is  given  by  Eq.  (1)  and  the 
theory  of  uncertainty  propagation  [18].  An  approxi- 
mation on  a Taylor  series  expression  gives  the  rela- 
tionship between  the  variables  as 

ST/T  ~ 8P/4P  -I-  8ri/2ri  + hr^ri  + 85/25  . (11) 


The  first  term  on  the  right  hand  side  of  Eq.  (11) 
shows  the  relationship  to  the  contribution  from  the 
uncertainty  in  the  radiant  power  which  has  both 
Type  A and  Type  B components  as  discussed  above. 
The  last  three  terms  in  Eq.  (11)  show  the  relation- 
ship to  the  contributions  from  the  uncertainties  in 
the  measurement  of  the  geometrical  quantities,  n, 
rz,  and  5.  A summary  of  the  percentage  uncertain- 
ties in  the  geometrical  measurements  (Table  2)  and 
their  propagated  uncertainties  for  radiance  temper- 
ature measurements  are  given  in  Table  6. 

The  predicted  NIST  radiance  temperature  from 
the  calibration  equation  for  each  average  value  of 
temperature  measured  by  PRT  sensor  1 is  given  in 
Table  7 for  the  BB  calibration.  The  uncertainty 
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Table  6.  Uncertainties  in  the  measurements  of  geometrical  quantities  for  the  BB 


Quantity 

Uncertainty  in 
measurement  (1  cr) 

Propagated  uncer. 
for  temperature 
measurement  (1  cr) 

1) 

Distance  between  apertures  (&S/5) 

0.136% 

0.068% 

2) 

Radius  of  ACR  aperture  (SrJrj) 

0.003% 

0.002% 

3) 

Radius  of  BB  aperture  (Sri/ri) 

0.2% 

0.1% 

TotaP  0.12% 

“ Square  root  of  sum  of  squares  of  1,  2,  and  3 given  above. 


Table  7,  Predicted  radiance  temperatures  and  uncertainties  for  BB 


Nominal 

Sensor  1 
data 

(K) 

Predicted 

radiance 

temp 

(K) 

ACR 

char. 

(K) 

Systematic  uncertainties 
Geometry  Diff.  Total 

meas.  cal.  b 

(K)  (K)  (K) 

1 cr  random 
uncertainties 

(K) 

Expanded 
uncer.  (U) 

(K) 

200 

199.906 

200.3 

0.06 

0.24 

0.07 

0.25 

0.10 

0.7  (0.4%) 

250 

224.777 

225.4 

0.07 

0.26 

0.07 

0.28 

0.08 

0.7  (0.3%) 

250 

249.678 

250.4 

0.08 

0.30 

0.08 

0.32 

0.07 

0.7  (0.3%) 

275 

274.683 

275.6 

0.09 

0.33 

0.05 

0.35 

0.05 

0.8  (0.3%) 

300 

299.547 

300.7 

0.09 

0.36 

0.06 

0.38 

0.04 

0.8  (0.3%) 

325 

324.429 

325.8 

0.10 

0.39 

0.07 

0.41 

0.04 

0.8  (0.3%) 

350 

349.362 

350.9 

0.11 

0.42 

0.07 

0.44 

0.04 

0.9  (0.3%) 

375 

374.120 

375.8 

0.11 

0.45 

0.08 

0.47 

0.05 

1.0  (0.3%) 

400 

399.071 

400.9 

0.12 

0.48 

0.08 

0.5 

0.07 

1.0  (0.3%) 

components  due  to  systematic  effects  are  listed 
separately  based  on  Eq.  (11):  1)  the  uncertainty  in 
the  power  5P/4P,  due  to  the  characterization  of 
the  ACR  as  an  absolute  detector  [1]  is  listed  in 
column  4;  2)  the  uncertainty  in  defining  the  geome- 
try as  given  in  Table  6 is  shown  in  column  5;  and  3) 
the  uncertainty  in  diffraction  calculations  due  to 
approximations  is  shown  in  column  6.  The  total  sys- 
tematic component,  b,  which  is  the  root  sum  of 
squares,  is  given  in  column  7.  The  Type  A compo- 
nent, s{Th),  of  the  uncertainty  is  obtained  from 
least-squares  analysis,  and  the  expanded  uncer- 
tainty, U,  as  shown  in  column  9,  is  obtained  by  ex- 
panding Eq.  (9)  as  follows  [19]: 

U = [(k  + l)  F{95;k  + l,n-k-l)y>^ 

[s\n)  + h^r-  (12) 


The  multiplying  factor  in  Eq.  (12)  for  n = 21,  k = l 
for  the  BB  experimental  data  is  2.6.  The  value 
shown  in  parenthesis  in  column  9 is  the  expanded 
uncertainty  given  as  a percentage  of  the  measured 
temperature. 

The  statistical  procedure  that  leads  to  Eq.  (12) 
did  not  take  into  account  possible  correlations  be- 
tween Type  B uncertainties  because  of  the  relation- 
ship between  variables  shown  in  Eq.  (11).  However, 
calculations  of  expanded  uncertainty,  U,  using  the 
complete  covariance  matrix  and  least  square  fitting 
using  appropriate  weights  have  been  carried  out 
and  the  results  are  found  to  be  the  same  as  given  by 
Eq.  (12).  The  calibration  constants  in  Eq.  (10)  also 
are  found  to  be  essentially  same  as  reported  earlier. 

Figure  3 shows  the  temperature  value  from  PRT 
sensor  1 plotted  on  the  A'-axis  and  the  measured  ra- 
diance temperatures  (Table  5)  and  corresponding 
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predicted  radiance  temperature  (Table  7)  on  the 
y-axis.  The  solid  line  in  Fig.  3 connects  the  pre- 
dicted radiance  temperature  values  and  represents 
Eq.  (10).  Figure  4 shows  the  95%  confidence  band 
for  the  difference  between  the  calibrated  tempera- 
ture and  the  PRT  temperature.  The  expanded  un- 
certainty given  in  Table  7 is  within  0.5%. 


Fig.  3.  Measured  temperatures  and  calibrated  radiance  tem- 
perature from  polynomial  fit  plotted  as  a function  of  blackbody 
temperature  setting  using  PRT  sensor  1. 


Fig.  4.  Solid  line:  Difference  between  the  deduced  radiance 
temperature  and  the  temperature  measured  by  the  blackbody 
PRT  sensor  1 plotted  as  a function  of  blackbody  temperature 
setting.  The  points  represent  the  same  for  measured  radiance 
temperatures  and  are  not  all  of  equal  weights.  Dashed 
lines:  Upper  and  lower  bounds  of  95%  confidence  bands. 

However,  it  is  always  an  open  question  why  for 
example,  at  the  nominal  blackbody  temperature  of 
the  400  K does  the  temperature  sensor  1 measure 
399.071  K whereas  the  radiance  temperature  is 


400.9  K as  shown  in  Table  7.  The  calibration  curve 
does  not  address  the  effects  of  systematic  errors 
such  as  unaccounted  scattered  light  entering  the 
ACR,  aperture  warming  of  the  blackbody  etc.  How- 
ever, experimental  designs  to  eliminate  such  sys- 
tematic effects  are  planned  and  implemented 
within  the  time  allotted  for  customer  calibrations. 
The  experimental  checks  show  these  effects  should 
be  negligible.  For  future  calibrations,  spectral  in- 
strumentation is  also  being  added  to  address  this 
question  further  by  measuring  the  emittance  from 
the  blackbody  as  a function  of  wavelength  between 
2 p,m  to  30  p,m  and  comparing  it  with  the  currently 
assumed  value  of  unity. 

4.  Summary 

The  LBIR  facility  is  in  operation  for  customer 
blackbody  calibrations  using  the  ACR  as  the  abso- 
lute detector.  The  experimental  procedure  and 
analysis  of  data  that  generate  the  calibration  report 
for  the  customer  are  described  in  detail,  and  a re- 
cent blackbody  calibration  at  the  facility  serves  as 
an  example  of  the  procedure.  Calibration  uncer- 
tainties in  radiance  temperatures  of  less  than  1% 
have  been  realized  for  typical  flux  levels  of  25  nW 
and  above  at  the  ACR  aperture.  Experimental  pro- 
cedures to  reduce  systematic  uncertainties  are  con- 
tinuously investigated  and  implemented  to  improve 
the  accuracy  of  measurements.  Future  improve- 
ments of  noise  control  in  the  ACR  operation  could 
allow  measurements  down  to  10  nW  at  the  ACR 
aperture  with  1%  uncertainty.  Capability  for  mea- 
suring spectral  emittance  of  customer  sources  in  the 
2 to  30  micrometer  wavelength  region  is  being  de- 
veloped. 
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Disclaimer 

References  made  in  this  paper  to  particular 
brand  names  or  specific  suppliers  of  a service  are 
made  for  ease  of  understanding  by  the  reader  and 
do  not  constitute  an  endorsement  of  products  or 
service  by  the  National  Institute  of  Standards  and 
Technology  over  other  competitive  suppliers  of 
similar  products  or  service. 
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